Optimal probabilistic estimation of quantum states 
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We extend the concept of probabilistic unambiguous discrimination of quantum states to quantum 
state estimation. We consider a scenario where the measurement device can output either an 
estimate of the unknown input state or an inconclusive result. We present a general method how 
to evaluate the maximum fidelity achievable by the probabilistic estimation strategy. We illustrate 
our method on two explicit examples: estimation of a qudit from a pair of conjugate qudits and 
phase covariant estimation of a qubit from A'' copies. We show that by allowing for inconclusive 
results it is possible to reach estimation fidelity higher than that achievable by the best deterministic 
estimation strategy. 
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I. INTRODUCTION 

The laws of quantum mechanics impose fundamental 
bounds on the amount of information that can be ex- 
tracted from the measurements on quantum states. In 
particular, it is not possible to exactly determine an un- 
known quantum state from a single copy. More gener- 
ally, nonorthogonal quantum states cannot be perfectly 
deterministically discriminated if only a single copy of 
the state is available. The optimal estimation and dis- 
crimination of quantum states has attracted a lot of at- 
tention during recent years. This interest was largely 
stimulated by the rapid development of quantum infor- 
mation theory. Quantum measurement forms an essen- 
tial ingredient of practically every protocol for quantum 
information transmission and processing since it converts 
the quantum information carried by a quantum system 
onto classical information. 

Two different strategies were proposed in the literature 
to optimally discriminate among nonorthogonal quantum 
states. The first option is to minimize the discrimination 
error, i.e. the probability of a wrong guess 
In this case the measurement device always provides a 
guess of the state, which sometimes may be wrong. An 
alternative approach pioneered by Ivanovic, Dieks and 
Peres 0, H is the unambiguous probabilistic discrim- 
ination which allows for perfect error-free identification 
of non-orthogonal states at the expense of a fraction of 
inconclusive results. It was shown that n pure quantum 
states from a set can be unambiguously dis- 

criminated if and only if they are linearly independent 
0, 0, ^> Q • Recently, the concept of unambiguous 
discrimination was extended to mixed quantum states 
^4, 15, 16, 17, 18]. More generally, it was shown that by 
allowing for a fraction of inconclusive results it is possible 
to reduce the probability of a wrong guess even if it is not 
pos sible to achieve the perfect error-free discrimination 

The optimal quantum state estimation can be thought 
of as a limiting case of quantum state discrimination 
among infinitely many states forming a continuous set. 
The canonical example of the state estimation task is a 



determination of an unknown state of a qubit from a sin- 
gle copy if a-priori the state could lie anywhere on the 
surface of the Bloch sphere The similarity of the 

(mixed) estimated state p^st with the true pure state \ip) 
is quantified by the fidelity F = (^/;|pest|V') ^'Hd the opti- 
mal state estimation strategy is defined as the one which 
maximizes the average estimation fidelity. During recent 
years optimal estimation strategies were established for 
a wide class of input sets of states including optimal uni- 
versal estimation of qubits [23, |^ and qudits , opti- 
mal phase-covariant estimation of qubits and qudits 
[2^ and optimal estimation of coherent states Also 
optimal estimation of mixed states 25] has been studied. 

In this paper we generalize the concept of unambiguous 
state discrimination to state estimation. We consider a 
scenario where the measuring apparatus can either out- 
put an estimate of the state or an inconclusive result. 
We shall show on explicit examples that with such prob- 
abilistic estimation strategy it is possible to increase the 
estimation fidelity above that achievable by the optimal 
deterministic estimation. 

The rest of the paper is structured as follows. In Sec- 
tion II we establish a general formalism for the determi- 
nation of the optimal probabilistic covariant estimation 
strategy. We shall show that the maximum attainable 
fidelity can be evaluated as a maximum eigenvalue of a 
certain operator. In Section III we apply the method to 
the determination of optimal probabilistic estimation of 
a qudit state \4>) from a pair of complex conjugate qudits 
This represents an extension of the well stud- 
ied problem of optimal estimation of a qubit from a pair 
of orthogonal qubits |2^ to d-dimensional quantum 
system. Recently, the optimal deterministic estimation of 
from a single copy of \4>)\4'*) was addressed by Zhou 
et al. who numerically calculated the maximum achiev- 
able fidelity We shall show that this numerically 
obtained fidelity is actually the fidelity of optimal prob- 
abilistic estimation and we will find a simple analytical 
formula for it. We will also derive from the first princi- 
ples the maximum deterministic estimation fidelity -Fdct 
for this case. Remarkably, -fdot turns out to coincide with 
the fidelity corresponding to the analytically found local 
extremum point in Ref. j28l ]. As a second example we 
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shall consider in Section IV the optimal phase covariant 
estimation of a qubit from N copies. We shall show that 
for N > 2 the probabilistic estimation strategy achieves 
strictly larger fidelity than the deterministic one and we 
shall compare the asymptotic behavior of the fidelities 
for large N. Finally, Section V contains conclusions and 
a brief summary of the main results. 



II. OPTIMAL PROBABILISTIC ESTIMATION 

Consider a set of pure quantum states j^* (?/>)) which 
are parametrized by a state {ip) and let S denote the set 
of all admissible \ip). The a-priori probability distribu- 
tion of lip) labeled by dip satisfies J^dip = 1. Using this 
notation we can treat in a unified way more complex situ- 
ations such as the estimation of \ip) from N copies of the 
state, when \'^) = \tp)'^^ . The goal of the quantum state 
estimation is to determine the state {ip) as precisely as 
possible by performing a generalized quantum measure- 
ment described by a positive operator valued measure 
Il{(p)d(p on Here \(p) e S and the detection of Il{(p) 
implies that the state \(p) is given as the estimate. The 
optimal measurement strategy generally depends on the 
set S and on the a-priori probability distribution dip. 

In the probabilistic estimation one allows for incon- 
clusive results where the machine does not produce any 
estimate of the state. This null outcome is associated 
with a POVM element IIq and the whole POVM should 
satisfy the completeness condition: 



/ n(</))d0 + no = 1, 



(1) 



where 1 is the identity operator on the Hilbert space 
spanned by the states {"^{ip)). The success of the esti- 
mation procedure can be conveniently quantified by the 
average fidelity. Consider first a particular input state 
IV'). The normalized fidelity of the estimation of 1-0) can 
be expressed as 



F(^) 



1 



{^{i;)\UicP)\'if{i;))\{m\'d<P, (2) 



where 



p(^) = / {^mii{^)\^mdcp (3) 

Js 

is the probability of producing an estimate of the state 
and, consequently, 1 — P{')p) is the probability of incon- 
clusive outcome. We choose as a figure of merit the nor- 
malized average fidelity. 



F = 



J^F{^j)P{^)d4, 
JsPWdi; 



(4) 



The maximization of is a complicated task since in 
general it requires the optimization of infinitely many 



POVM elements 11(0). The problem simplifies consid- 
erably if the states \ip) and |^ ('(/')) form orbits of some 
group G, S = G, and if dip is an invariant measure in- 
duced by the Haar measure on G. In the rest of the 
paper we will assume that this is the case. We then have 
10) = C/(0)|O) and |*(0)) = V^(0)|*(O)), where C/(0) 
and Viip) denote unitary representations of the group G. 
With slight abuse of notation we use tp to label the ele- 
ments of G. It can be shown that due to the underlying 
group structure the optimal POVM which maximizes F 
can always be chosen to be covariant and all the POVM 
elements are generated from a single element. 



Uc{cp)^v{(P)ncV\^), 



(5) 



and He = nc(0). 

On inserting the expression Q into Eq. ^ we obtain 



F = 



Tr[^nc] 



where 



and 



|*(7^))(vE'(0)||(0|O)|2d0 



I |vl/(0))(vl/(0)|d0. 
Ig 



(6) 



(7) 



(8) 



Note that the expression ^ is formally similar to the 
formula for the fidelity of the optimal probabilistic com- 
pletely positive map that approximates some unphysical 
operation which was derived in Ref. Upon intro- 

ducing 



fie = a'/^UcA'/^ 

we can rewrite Eq. © as 



F 



Tr[A-i/2i^^-i/2n^] 



Tr[n 



(9) 



(10) 



c\ 



It follows that the fidelity is bounded from above by 
the maximum eigenvalue /Xmax of the operator M = 
A-^/^RA-^/^. If the Hilbert space spanned by |^(0)) 
is finite dimensional then there exists a POVM which at- 
tains the maximum fidelity -Fmax = /^max and produces 
the estimate of |0) with nonzero probability P > 0. Let 
l^max^ , j = 1, . . . , J, be the eigenvectors corresponding to 
the maximum eigenvalue yitmax- 

Then the POVM element 
Ilcopt which generates the optimal covariant POVM can 
be expressed as 

,7 

Hcopt = ^ ^,fe|^--)(^--|^-l/2. (11) 



The coefficients 7r., i. must be chosen such that 11, 



n, 



t 

C,opt 



Ilcopt and 



C.opt 



Vii:)Uc,optV\iP)d4> < 1. 



> 0, 



(12) 
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The coefficients iTjk may be optimized such as to maxi- 
mize the average probabihty of success P under the con- 
straints Ilcopt > and (|12|l . This is an instance of 
a semidefinite program which is a convex optimization 
problem that can be very efficiently solved numerically 

If the maximum eigenvalue is non-degenerate, J = 1, 
then the optimal POVM element Hcopt is proportional 
to a rank-one projector 



n 



C,opt 



-1/2 



(13) 



and the normalization constant Af has to be chosen such 
that (O holds. 



III. PAIR OF CONJUGATE QUDITS 

A. Optimal probabilistic estimation 

In this section we will investigate the optimal proba- 
bilistic estimation of a state of a single qudit j?/') from 
a pair of conjugate qudits, |^'(?/')) — IV')IV'*)- The a- 
priori distribution d-ip is assumed to be induced by the 
Haar measure on the group SU(d). This scenario is a 
generalization of the estimation of a qubit from a pair 
of orthogonal qubits .27] to dimensions d > 2. The 
corresponding operators A and R can be easily evalu- 
ated with the help of the Schur Lemma. The unitary 
representation JJ^^ of the group SU((i) acts irreducibly 
on the totally symmetric subspace of N qudits. Taking 
into account that \ip*){ip*\ = {\ip){ip\)'^ and exchanging 
the order of integration and (partial) transposition we 
obtain 



A 



1 



d{d+l) 



(l + rf$+), 



where $+ = |$+)(<i>+| and 



i$+) = ^Eij>ij) 



(14) 



(15) 



is a maximally entangled state of two qudits. We also 
find that 



R 



1 



D+{d) 



Tr3[li® 12® 10)3(01 (n+3) 



T21 



(16) 



Here 11^23 is the projector onto the symmetric subspace of 
three qudits, D^{d) — ^d{d + l){d + 2) is the dimension 
of this subspace, T2 denotes partial transposition with 
respect to the second qudit, and Tv^ stands for the partial 
trace over the third qudit. After some algebra we find 
that R can be expressed as 



R 



1 



li 



d{d+l){d + 2) 
+Vd\00){^+\ + Hi (g) |0)2( 



\/rf|$+>(00| 

K|0>l(0|(g)l2 



(17) 



Instead of calculating the eigenvalues of A~^/^RA~^^^ 
we can equivalently look for the eigenvalues of RA~^ be- 
cause the eigenvalues of these two operators coincide and 
the latter is easier to deal with. The operator A can be 
easily inverted. 



= d{d+l) 



d+1 



(18) 



and we arrive at 
1 



RA-' = 



d+2 



I12 - 



-$+ + \/d|$+)(00| 



d+1 



d+1 

|00>($+| + li® 10)2(01 -I- |0)i(0|®l2 

(19) 



This operator possesses only four different eigenvalues 
which can be expressed analytically for arbitrary d. The 
eigenvalue /ii = l/((i -I- 2) is d{d — 2)-fold degenerate 
and the eigenstates read |j)i|fc)2, where j ^ 0, k ^ 0. 
The second eigenvalue ^4 — 2/{d + 2) is {2d — 2)-fold 
degenerate with eigenstates |0)i|j)2 and |j)|i0)2, j 7^ 0. 
Finally, the last two eigenvalues are non-degenerate and 
can be expressed as 



A*3,4 



d + 2 



1 ± 



2id+l) 



(20) 



The maximum eigenvalue is 113 for all d > 2 and the 
fidelity of the optimal probabilistic estimation of \'>p) form 
a pair of conjugate qudits is equal to this eigenvalue, 



max,prob 



d + 2 




2id+l) 



(21) 



The numerical values for the optimal estimation fidelity 
F± obtained by Zhou et al. [28|| fully agree with the above 
analytical formula so their global optimization actually 
yielded the optimal probabilistic estimation strategy. 

The optimal probabilistic covariant POVM is gener- 
ated by the POVM element Ilc.opt = kcopt) (7rc,opt|, 
where 



kcopt) « |oo) 



2d y/2{d+l) - Vd 



d+1 



d + 2 



(22) 



B. Optimal deterministic estimation strategy 

For d > 2 the optimal estimation strategy obtained 
above cannot be made deterministic and there is a 
nonzero probability of inconclusive results. Thus a ques- 
tion arises what is the optimal deterministic strategy of 
estimation of \ip) from a single copy of the state {ipip*}- 
When seeking an answer to this question we can restrict 
ourselves to the covariant POVMs. Since the probability 
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of inconclusive results should vanish, we have Hq = and 
the completeness condition for the POVM becomes 



U{i;) ® f7*(V')ncC/^(V') ® U^ii^)di; = 1. (23) 



Recall that U{ip) is a unitary acting on the Hilbert space 
of a single qudit and U{iIj)\0) — lip). We should max- 
imize the estimation fidelity F — Tt[IIcR] under the 
above completeness condition. In the present case this 
is equivalent to maximizing F under simpler constraints 
that can be obtained from H23() . In particular, by calcu- 
lating the trace of Eq. (|23|l and by taking into account 
the invariance U (g) C/*|<I'+) = 1$+) we find that 



Tr[nc$+] = 1 



(24) 



must hold. 

The constraints H24() can be accounted for by introduc- 
ing two Lagrange multipliers Ai and A2 and our task is 
to maximize 

J^lUc] = Tt[RI1c] ~ AiTrpc] - hTi'l^+Uc] (25) 

under the constraints (|24|l and Il r- > 0. This is an in- 
stance of a semidefinite program jSfl- For this class of 
convex optimization problems one can straightforwardly 
derive the extremal equation for the optimal He and we 
obtain 

(i?- All - A2$+)nc = 0. (26) 

Moreover, we also find the optimality condition, 

Ail + A2<I'+ -i?> 0. (27) 

If (|26|l and 1)27(1 hold simultaneously, then He is the opti- 
mal one which maximizes the fidelity. To prove this state- 
ment we take the trace of a product of Eq. H27() with an 
arbitrary lie which satisfies all the constraints imposed 
on it. We get Trpci?] < \id^ + A2 hence the Lagrange 
multipliers provide an upper bound on the achievable fi- 
delity which is saturated if the POVM satisfies Eq. (|26|l . 

The optimal POVM has qualitatively similar struc- 
ture as the optimal probabilistic POVM H22I) . namely, 
Iladct = |7i'c,dct)(7i'c,dot| where 

\ncAct)^Vd{d + T)m-{Vd + T^l)\<i>+). (28) 

By construction, this POVM satisfies the completeness 
condition . On inserting llc^dct h^to Eq. we can 
solve for the Lagrange multipliers and we get 



Ai 
A2 



did+l){d + 2) 



r^3 



d{d+l){d + 2) 



2d^ -2d -A 



dVd + T 



This choice guarantees that Eq. (|26|l holds for any d. 
To prove the optimality of the POVM (|28|l it remains to 
check the inequality Since Ai > 2/[d{d+ l)(d + 2)] 

the only nontrivial part is the verification of the positive 
semidefiniteness of the operator in the two-dimensional 
subspace spanned by |00), Let 112 denote the pro- 

jector onto this subspace and consider the 2x2 matrix 
K = n2(Ai]l + A2^'+ - -R)n2. One eigenvalue of K is 
zero due to the optimality condition (|26|l . To prove that 
K > it thus suffices to show that TrK > and after 
some algebra we find 



TtK = 2Ai + A2 



d + 6 



d{d+l){d + 2) 



(30) 



It can be shown that TrK is a growing function of d and 
that it is positive for all integer d > 2. This concludes 
the optimality proof. 

The fidelity of the optimal deterministic estimation 
corresponding to the optimal covariant POVM reads 



Fr 



■dot 



d^{d + 2) 



Sd^ - 4d + 4 



2d^ + 2d - 4 



Vd + T 



(31) 

This expression agrees with the formula given by Zhou et 
al. [23 • In that paper, the authors claimed that this is 
only a local maximum of the fidelity and they calculated 
the global maximum of the fidelity numerically. Our find- 
ings provide a precise interpretation of their results. The 
local maximum is in fact the maximum achievable fidelity 
of deterministic estimation from a pair of co njug ate qu- 
dits while the global maximum given in Ref. |28| | corre- 
sponds to the optimal probabilistic estimation strategy 
which allows for inconclusive results. 

It is instructive to explicitly evaluate the probability 
that the machine outputs an estimate |0) for an input 
state \tpip*). We have 



P^(0|^) = i 



dVdTT|(V'|o)p -Vd + T+i 



(32) 



(29) 



Note that this probability is zero if the overlap of the 
true state j?/;) with the estimated state |0) is equal to 
\{iIj\0)\^ = {Vd+1 - l)/(d\/c?+ 1). It is interesting to 
compare this with the optimal estimation from a pair 
of identical qudits IV'V')! where the optimal covariant 
POVM is generated by |7rc,||) = ^^(d -t- 1)/2|00) and 
the corresponding probability of guessing |0) for the in- 
put state \'ip) reads P\\{QW) = \d(d + which 
vanishes only if the state j-i/') is orthogonal to |0). In par- 
ticular, for d = 2 the probability P\\ = only if \tjj) — |1) 
while P± is zero for all states on a certain circle of the 
Bloch sphere. This observation gives some more insight 
into why the state {ip) can be estimated with higher pre- 
cision from IV'V'*) than from l^ijj). 
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IV. OPTIMAL PROBABILISTIC ESTIMATION 
OF EQUATORIAL QUBITS 

In this section we will investigate the optimal phase- 
covariant probabilistic estimation of a qubit. We shall 
assume that it is a-priori known that the qubit state is 
located on the equator of the Poincare sphere, \ilj{(p)) ~ 
-^(|0) + e^'^\l)). The state is thus characterized by a 
single parameter - the relative phase (p. Starting from the 
seed state :^(|0) + |1)) all states \fp{(p)) can be obtained 
as the orbit of the Abelian group U(l) which generates 
rotations of the Bloch sphere about z axis. To make our 
treatment general we will consider optimal estimation of 
liplip)) from N input copies. 

This scenario corresponds to a typical phase-shift mea- 
surement, where N particles pass through an interferom- 
eter which applies an unknown relative phase shift ip to 
one of the states of the particle and the goal is to de- 
termine if as precisely as possible. The optimal phase- 
estimation strategies which reach the so-called Heisen- 
berg limit Ag9 « require entangled input states of 

particles [3ll. Here we show that even for the prod- 
uct state I'tpif))'^^ it is possible to probabihstically im- 
prove the precision of ip estimation, as witnessed by 
the improved asymptotic scaling of the optimal fidelity 
oh cc Tri- However, it should be noted that 
this apparent improvement is achieved only for the sub- 
ensemble of conclusive measurement outcomes while the 
inconclusive outcomes are neglected. 

The input state belongs to the N+1 dimen- 

sional symmetric (bosonic) subspace of the Hilbert space 
of N qubits and it can be written as follows, 




ikip 



\N;k). (33) 



Here \N]k) denotes a normalized fully symmetric state 
of N qubits with k qubits in state |1) and N — k qubits 
in state |0). 

The operators A and R can be determined from the 
formulas ||SJ| and (TJ, where we have to integrate over the 

phase shift (p, J d^jj — Jq^ ^'^'■P- After the integration 
we obtain 



A 



N 



1 ^ 

fc=0 



\N-k){N-k\ 



(34) 



0.01 



0.008 



0.006 



AF 



0.004 



0.002 




FIG. 1: The difference of the fidelities of optimal probabilistic 
and deterministic estimation /S.F — i^max.prob — finax,dct is 
plotted as a function of the number of copies A*' of the qubit 
state \tp{ip)). 



imum eigenvalue determines the maximum achievable es- 
timation fidelity can be easily calculated and we have 

M = A-^I^RA-^I^ = + i Y.^^^ + 4)- (36) 

Instead of directly working with M let us consider 
the operator M — AM — 21 and the eigenvalues 
of M are then related to the eigenvalues p,j of M by 
/ij = {jlj + 2)/4. The matrix M is tridiagonal and 
its characteristic polynomial is given by the Tcheby- 
chev polynomial of the second kind, det(M — Al) 
UN+ii^^)- Maximum eigenvalue of M is thus given by 
the largest root of Un and with the help of the definition 
t/jv(cos6') = sm[{N + l)6']/sin6l wc arrive at [H 



— 2 cos 



iV + 2 



(37) 



The maximum fidelity can be determined as (/imax + 2)/4 
and we finally obtain 



1 



max.prob 



1 + cos 



N- 



(38) 



and 



It is instructive to compare the fidelity i^max.prob with 
the fidelity of the optimal deterministic phase covariant 
estimation of a qubit from N copies 21], 



R = 



1 ^ 

N+l X! 



N 



2N+1 ^ y 



\N-k)(N-k\ 




k=l 



(Xfc + 4), (35) 



where Xk — \N; k) {N; k — 1|. Since the operator A is 
diagonal in the basis |A^; fc), the operator AI whose max- 



;,dot 




k=l 



(39) 



It follows that for = 1 and = 2 Fmax,prob = i^max.dct 

hence it is not possible to improve the fidelity of estima- 
tion by allowing for some fraction of inconclusive results. 
However, if A^ > 3 then Fmax.prob > ^Inax.det and the op- 
timal probabilistic estimation strategy attains a strictly 
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FIG. 2: Probability of successful estimation P versus the 
number of copies TV of the qubit state \tp{(p)). 



larger fidelity than the optimal deterministic strategy. 
This is illustrated in Fig. 1 which shows the difference 
between the fideHties and 

Since the maximum eigenvalue of M is non-degenerate, 
the optimal covariant POVM is generated by the rank- 
one projector The maximization of the success 
probability P = Tr[nc^opt^] is equivalent to the mini- 
mization of J\f under the constraint (|12l) . Let |/Umax) = 
Sfc=o '^k\N; k) be the normalized eigenvector of M with 
the eigenvalue /imax- Then it is optimal to choose 



A/" = 2^ max 



(40) 



Numerical calculation reveals that P decreases exponen- 
tially with growing TV, see Fig. 2. 

The further facilitate the comparison of the fidelities 
let us analyze their asymptotic behavior for large TV. For 
the probabilistic estimation strategy we obtain 



max,prob 



1 - 



1 



4 (TV 



(41) 



and we can write 1 — Fi„ax,prob = 0{N~'^). On the 
other hand, for the deterministic estimation we find that 
1 — -Fi„ax,dct ~ 0(TV~^). We can see that with growing TV 
-Finax.prob convcrgcs to Unity much faster than Fmax.dot- 
This superior scaling is achieved at the expense of a de- 
creasing probability of successful estimation P. 



V. CONCLUSIONS 

In the present paper we have generalized the concept of 
unambiguous quantum state discrimination to the realm 
of quantum state estimation. We have shown that by 
allowing for inconclusive results it is possible to increase 
the fidelity of estimation evaluated for the sub-ensemble 
of conclusive outcomes of the estimation process. We 
have established a general formula for the maximum fi- 
delity achievable by the probabilistic state estimation. 
The method was illustrated on two explicit examples. 
First, we have studied the optimal estimation of a qudit 
from a pair of conjugate qudits and we have provided an 
exact interpretation of the results recently obtained by 
Zhou et al. |2^. As a second example we have investi- 
gated the phase covariant estimation of a qubit from TV 
copies of the state. The present quantum-state estima- 
tion scheme could find applications in quantum commu- 
nication and it may potentially help to probabilistically 
improve the sensitivity and precision of measurements 
performed at the quantum limit. 
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